Van der Waals-Casimir-Polder interaction of an atom with a composite surface 
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We study the dispersion interaction of the van der Waals and Casimir-Polder (vdW-CP) type between a neutral 
atom and the surface of a metal by allowing for nonlocal electrodynamics, i.e. electron diffusion. We consider 
two models: (i) bulk diffusion, and (ii) diffusion in a surface charge layer. In both cases the transition to a 
semiconductor is continuous as a function of the conductivity, unlike the case of a local model. The relevant 
parameter is the electric screening length and depends on the carrier diffusion constant. We find that for distances 
comparable to the screening length, vdW-CP data can distinguish between bulk and surface diffusion, hence it 
can be a sensitive probe for surface states. 

PACS numbers: 34.35.+a - interactions of atoms with surfaces; 31.30.jh - long-range QED interactions. 



I. INTRODUCTION 

Recent progress in the understanding of the van der Waals- 
Casimir-Polder (vdW-CP) force between an atom and a sur- 
face allows by now to distinguish surface properties with re- 
spect to charge transport. Data on the temperature and atom- 
surface distance dependence provide excellent tools for such 
analysis. In particular the experiments on fused silica |ll have 
demonstrated a temperature dependence of the vdW-CP inter- 
action. Fused silica is considered as a dilutely doped semi- 
conductor that has a finite conductivity a. The data was fitted 
successfully to the potential for a dielectric surface which dif- 
fers from the one for a perfectly reflecting mirror, as consid - 
ered by Casimir and Polder J2l ■ In fact, since u/lo diverges at 
zero frequency, any nonzero conductivity will reproduce the 
perfect reflector result fl. It has been suggested to resolve 
this puzzle with the help of a nonlocal model of the elec- 
tromagnetic response, where charge diffusion and screening 
become essential at low conductances Q], although the ap- 
proach was met with criticism J5|. In a related experiment on 
the (macroscopic) Casimir force between a gold-coated sphere 
and a single-crystal silicon membrane fl, the results are con- 
sistent with a dielectric behavior in its pristine form (ignoring 
the a/oj tail). The significant change in the charge carrier 
density after laser illumination leads to a metallic response in 
local form, again in agreement with Casimir force measure- 
ments. Nonlocal theories have been worked out to understand 
the crossover between these limits flHUt, although experi- 
mental data favor a local description where the contribution of 
free charge carriers is omitted in the dielectric state . For 
the discussion whether nonlocal electrodynamics may be ap- 
plied to macroscopic Casimir interactions, and its consistency 
with thermodynamics, see RefsJ5ll9l 4l ill 1311 . 

The ability of probing charge transport may provide for 
an increased understanding of surface science by using the 
vdW-CP interaction as a probe. This viewpoint may be traced 
back to the seminal paper by Zaremba and Kohn jl-fl where 
the van der Waals potential (neglecting retardation) was cal- 
culated with a microscopic description of the many-electron 



response of a metal. Their analysis yields an expression for 
the reference plane with respect to which the atom-surface 
distance z is actually calculated. See Ref. fl5ll for a review 
of related methods. Dorofeyev IU6I1 analyzed the van der 
Waals (non-retarded) regime with the help of a nonlocal (k- 
dependent) electromagnetic response based on the surface 
impedance work by Kliewer and Fuchs (l7\\ ; the electrons 
were assumed to reflect specularly from the inner surface. 

In the present work, we compare two non-local models that 
can be understood as mesoscopic extensions of the work by 
Zaremba and Kohn: the first one allows for bulk diffusion 
as in Refs.J3,0], we denote it "continuous charge" (CC). In 
the non-retarded limit, this reduces to the Kliewer and Fuchs 
approach for a hydrodynamic dielectric function in the bulk. 
The model (ii) allows for diffusion only within a surface layer 
("charge layer" or CL), the bulk charges responding with a 
local conductivity. Such composite surfaces with charges in 
the bulk and on the surface, are fairly common in metallic 
systems. These surfaces are either covered with adsorbates 
or nanostructures and can, ej* be used as sensitive chemi- 
cal sensors and biosensors 11811 . or are disordered with quan- 
tum well states at the surface lfl9ll . Further motivation for a 
two-type charge model comes from studies of the anomalous 
heating of cold ions observed in miniaturized Paul traps that 
invoke surface charge fluctuations on the metallic electrodes 
1I20I- 24 1 . Composite surfaces have also been explored regard- 



ing surface plasmons and lead to a wide range of dispersion 
relations, as observed in different systems, see Ref. n2.5h and 
references therein. 

The hallmark of the non-local theory is charge diffusion. In 
the CC model, it is described by the diffusion coefficient D in 
the bulk, while D s describes surface diffusion in the CL. We 
identify a screening length uq for CC (cgs units) 
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where er, 



are the bulk conductivity and background di- 
electric constant, respectively (replace D with D s in the CL 
model). As the atom-surface distance z becomes compara- 
ble to the length scale ap, the vdW-CP interaction changes its 
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behaviour. We consider different materials where ao can be 
compared with the other two important length scales of the 
vdW-CP potential: the radiation wavelength Xa = c/f2, 
being the atomic transition frequency, and the thermal photon 
wavelength At = ac/fe^T at temperature T. As common for 
electric dipole transitions in the visible range, we assume here 
Xa -C At, but explore otherwise the full range of distances. 
After defining the charge models (section|II]i, we study in sec- 
tion|III]the case that ao is in the van der Waals (short distance) 
regime, ao <C Xa, focusing on the T = case. In section HV1 
the length ao is in the Casimir-Polder (intermediate) regime, 
A A <C ao -C At, while in section [V] we consider long dis- 
tances ao ^> At where screening affects the Lifshitz (ther- 
mal) regime of the interaction. In all these cases we find that 
the crossover from the vdW-CP potential for a local conductor 
to a dielectric occurs at z « ao. In addition, the dispersion in- 
teraction is a good surface probe in the sense that the crossover 
is different in the CC and CL systems. This difference is par- 
ticularly visible when the two limits, metallic and dielectric, 
are well separated. The relevant parameters are discussed in 
the conclusion. 



H. MODEL 

A. Atom-surface interaction potentials 

We use in this paper the general formulation of Wylie and 
Sipe ll26ll for the Casimir-Polder potential of an atom with a 
surface. Assuming the surface and the electromagnetic field 
in equilibrium at temperature T, the free energy of interaction 
is given by 



(2) 



which is a sum over Matsubara frequencies = 2-Klk B T /% 
along the imaginary frequency axis. The primed sum is tak- 
ing the I = term with a factor |. The atomic polarizabil- 
ity a.ij (to) is given in Eq.© below, the retarded Green tensor 



Gij (r, r; uj) is made explicit in Eqs. (1617b . In the limit T — > 
F reduces to the interaction energy 



U{v) = 



<li"V<>,/: /XX •'.„•:>. r; i£) 



2ir 



(3) 



This formulation applies for an atom in the ground state, for 
which the polarizability tensor is given by 
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where Vl eg is the transition frequency between the ground state 
(<?) and an excited state (e) with an electric dipole matrix el- 
ements g <H> e and df e = (g\di\e). We focus on a single 



resonance and assume rotational symmetry, so that tt e 



n 



and the polarizability is isotropic, a. 



a Si 



The electromagnetic Green tensor Gjj(r,r';w) in Eq.© 
provides the electric field Ei (r) radiated by a test dipole lo- 
cated at r' and oscillating with amplitude dj at frequency u>: 



(5) 



For a source outside a polarizable body, this electric field can 
be calculated within macroscopic electrodynamics (see, e.g., 
Ref ill ) and involves the reflection (or scattering) amplitudes 
of the body. It turns out that these amplitudes are sufficient to 
determine the Casimir-Polder interaction. The subtraction of 
the free-space part of the Green tensor (Lamb shift) is under- 
stood in the following. At a planar surface, only two principal 
polarizations p = TM, TE are relevant, and the reflection am- 
plitudes depend on frequency and a wave vector h parallel to 
the surface r p = r p (u>, k). As we put r = r' in Eqs.(f2] |3j, 
the planar symmetry implies that the Green tensor is diagonal 
with elements ll26Tl (cgs units) 
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the element G yy being identical to G xx . The vacuum decay 
constant for a field mode with frequency oj and parallel wave 
vector fc is 
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= v/fc 2 - L0 2 /C 2 
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with the root chosen such that Rc vq > and Im vq < 0. The 
reflection amplitudes are collected in Table |I]for the different 
surface models considered in this paper. The general Eq.((2), 
valid for any T, now takes the form 
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Let us briefly recall the assumptions behind the Wylie and 
Sipe approach 126J: the interaction energy is calculated in per- 
turbation theory with respect to the atom-field coupling, start- 
ing from a well-defined atomic level (here, the ground state). 
The temperature provides Boltzmann weights for the excited 
states of the electromagnetic field, thus including the interac- 
tion with blackbody radiation and its modification by the sur- 
face. The thermal population of excited states of the atom 
is negligible provided the Bohr frequency is large enough, 
fifl ^> k B T. Otherwise, a temperature-dependent polariz- 
ability should be used in Eqs.(f2] [9}. Finally, the surface re- 
sponse is worked out ignoring the presence of the atom and 
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assuming a linear response of the surface to electromagnetic 
radiation, consistent with common practice in surface spec- 
troscopy. By inspection of Eq.©, one notes that wave vec- 
tors up to k ~ 1/z are relevant for the interaction potential. 
At distances z much larger than the size of the unit cell, a 
macroscopic treatment of the surface response is therefore jus- 
tified. There have been discussions what kind of electromag- 
netic response may be used consistently for atom-surface po- 
tentials and macroscopic Casimir interactions in general, see 
Refs.j5l|^- n, 131. The approach of Ref.1261 is based on the 
fluctuation-dissipation theorem 027I1 . 



B. Surface response with charge diffusion 

Surfaces covered with thin layers of strongly localized 
charges or adsorbates have been studied in much detail in sur- 
face physics. For a general theory of their electromagnetic 
response, see Refs. 11281 12911 . We consider here a model intro- 
duced in Ref . 12211 where the surface is covered by a charge 
sheet with a charge density 7(2, y) (localized in the plane 
z = 0). The details of the electromagnetic response are 
worked out in Appendix lAl The sheet current in the layer re- 
sponds by diffusion 



J = -D,(w)Vi|7 
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where D s (us) = D a j (1 — ilot s ) is the surface diffusion coef- 
ficient and t s a surface relaxation time. The gradient appear- 
ing here is parallel to the layer. A surface conductivity term 
proportional to E is neglected, as justified for a small layer 
thickness, see Appendix IA 31 Charge conservation yields 
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so that the bulk current just below the layer, j z (0~), provides 
the influx into the surface layer. We take the bulk current re- 
sponse in the usual Ohmic form 
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where the Drude conductivity is a(uj) = u/(l — iuir) with a 
scattering time r. Eqs. d 1 01 1 21 > define the CL (charge layer) 
model. 

The CC (continuous charge) model, in contrast, is defined 
by the bulk charge density that has a diffusion constant D, i.e. 
Eq. ([T2I 1 is replaced by 



j = cr(w)E - D(w)Vp 
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The ("additional") boundary condition for the current is then 
jz(0 ) = 0, since there are no surface charges. 

The resulting reflection amplitudes are summarized in Ta- 
ble H] and more details on their derivation are given in Ap- 
pendix A. One notes that the TE polarization is not affected 
by the composite structure of the surface. This can be under- 
stood from the fact that surface charges are created by electric 



fields perpendicular to the surface, which are absent in this 
polarization. 

To illustrate the impact of the diffusive layer, we have cal- 
culated the local photonic mode density, i.e., the imaginary 
part of Gu(r, r;uj). This quantity can be measured from the 
spontaneous decay rate of an excited atom placed at r or 



from the heating rate of an ion trapped near a surface fl20 



2411 . The results shown in Figfl]illustrate the enhancement of 
the mode density at low frequencies (below the characteris- 
tic scale D( s )/z 2 left column). At large wave vectors (short 
distances, right column), there is a competition between ad- 
ditional modes (enhancing the mode density, CL model) and 
screening (reducing it, CC model). Note that for the parame- 
ters considered here the screening length do is much smaller 
than the diffusion length ~ ^D( s )/ui). An excited atom de- 
cays faster because diffusion along the surface broadens the 
field spot it creates, increasing the effective area where ab- 
sorption takes place. Calculations of the Casimir (plate-plate) 
interaction between materials with a nonlocal electromagnetic 
response have revealed qualitatively similar trends (compare 
Refs. SHI to Ref.JH). The experimental data of Ref.Jl 
are better described with a local rather than nonlocal theory, 
however, see Ref. [11] and the discussion in Refs. ll2lll3 l. 



III. VAN DER WAALS (NONRETARDED) REGIME 

This regime corresponds to short distances where retarda- 
tion is negligible 



z < X A = c/fi 



(14) 



with a typical value « 100 nm for transitions in the vis- 
ible range. The Van der Waals interaction follows a power 
law T(z) ~ 1/z 3 for a material with a local response (Drude 
metal or dielectric). We consider in this section the situation 
that the screening length (see Introduction) satisfies ao < Aa; 
this corresponds to electron densities typical for metals. 

We start from the zero-temperature expression for the inter- 
action potential, combining Eqs. (01!] [6] [7}: 
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where now Vq = k 2 + £ 2 /c 2 . The dominant ranges of the 
integrals are around k~l/z^>Cl/c, due to the exponential, 
and £ < ft, due to the polarizability a(i£). This allows to 
simplify Eq. (ITSb by taking £ <C cfc and v ~ k, so that for 
the CL model (see Table |I| 
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hydrodynamic bulk charge 


(CC) charge layer (CL) 
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vo + v 


wo — V 
vo + v 
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evo — v 
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evo+v + ie-e^)^ 


e a v + v 



TABLE I: The reflection coefficients vte(u, k) and ttm^, k) for two models for a nonlocal surface response. 

Symbols used: The dielectric function e(ui) is defined at dl8b while e s (ui,k) is defined at l |17| >. Spatial decay rates are for the vacuum 
vo(cv,k), Eq.©, for transverse fields in the medium v(io,k) — \Jk 2 ~ e(uj)uj 2 /c 2 with Re« > and for compressional charge waves 
vi(lu, k) — yjk 2 — io;e(a;)/[e 00 D(ii))] with Re vi > and D(uj) — D/(l — iujr). 




FIG. 1: (left column) Local photonic mode density near a metallic surface, described as a local Ohmic conductor (solid blue curves), by a 
continuous bulk charge (CC) model (dotted yellow curves, superimposed), and a diffusive charge layer (CL) model (dashed purple curves). 
Thin green lines: free space mode density. We plot the imaginary part of the Green functions in Eqs.(|6]|7}; top row: polarization parallel to 
surface (G xx ), bottom row: perpendicular polarization (G zz ). Parameters: conductivity a = 3.6 ■ 10 17 rad/s, typical for Al, relaxation time 
t = t s — 10~ 15 s, diffusion constant D s = D = 5 • 10 3 cm 2 /s, background permittivity = 1, giving a screening length ao ~ 0.3 nm, 
and distance z = 10 /im. The plots include the free space response that provides the dominant ~ u> 3 scaling at high frequencies. Characteristic 
diffusion frequency for these parameters: loo = D/z 2 — 5 ■ 10 9 rad/s (left arrow); right arrow: inverse photon round trip time lo c = c/ (2z). 
(right column) Distance dependence of the local mode density near a metallic surface, described by the same models as in the left column. 
Same parameters, except that the frequency is fixed to to = 10 rad/s (near-infrared). The plots include the free space response, that leads 
to a constant limit |cj 3 at large distance (thin horizontal lines). The oscillations are due to partial standing waves formed above the surface. 
Diffusion length for these parameters: Id = ^/2D/uj ~ 31 nm (left arrow); right arrow: reduced wavelength A = c/lj. 



The reflection amplitude involves a "surface dielectric func- 
tion" given by (see Ref. 12211 and AppendixlAl 



e s (uj, k) = Eoo + 
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where D s (lo) is the surface diffusion coefficient [see Eq.jTOb I. 
The conventional Drude dielectric function (local), with a bulk 



relaxation time r and a high-frequency asymptote e^, is 

. . 471W 

e(w) = £ 



uj(1 — \ujt) 



(18) 



which happens to be the limiting form of Eq.dTTIi when D s — > 
0. From the fc-dependence in e s (w,fc), we identify the di- 
mensionless ratio z 2 47rcr/(e 00 L' s ) = (z/ao) 2 that defines the 
screening length ao consistent with the estimate (fTJ of the In- 
troduction. For fcap <C 1, the reflection amplitude in Eq.JToTi 
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recovers the local behaviour, while at very short distances, 
kciQ ^> 1, the diffusive (nonlocal) term dominates. The latter 
case implies that the conductivity contribution is suppressed in 
Eq.dT7b. leaving only the background dielectric constant e^. 

If £oo > 1. both limiting cases (conductor and dielectric) 
show a van der Waals interaction that follows the familiar 
C3/2 3 power law, but with different C3 coefficients. For the 
dielectric, 

local did., z« Ax: ^(*) « ~ ^5" !°° 7 ? (19) 
while the local Drude conductor gives (2 



local metal, zCi 



U(z) 
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Here, uj 2 = 2na/r is the surface plasmon frequency, and we 
have neglected Ohmic losses (i.e., u 3 , ft 1/t). The com- 
posite surface models with their nonlocal response give a van 
der Waals interaction that crosses over between these two lim- 
its (see Figf5]below). This is similar to what has been analyzed 
at large distances by Pitaevskii 01 . 

We focus here on the somewhat academic case of a simple 
free-electron metal (coc = 1 ), where the nonlocal surface re- 
ponse changes even the exponent at short distances because 
the C3 coefficient in Eq.fO vanishes. To get insight into the 
small-distance behaviour, we expand the CL reflection coeffi- 
cient ttm at high momentum and get for r s = r 



k ~ l/o > C/c : 
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e s (i£, k)k — v 
e s (i^k)k + v ~ 1 + 2(fca ) 5 
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Note that the imaginary frequency £ drops out in this case and 
the suppression ~ l/k 2 on scales shorter than the screening 
length ao- The £ integral in Eq.dToll can then be performed, 
and one gets for z -C ao the simple result 



CL, z < a : U{z) 
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instead of the 1/z 3 power law. In Figj2] we compare the exact 
evaluation of the van der Waals potential (T5[ (dashed curve) 
to the numerical integration over k of the approximate re- 
flection amplitude (fJTJ (solid gray curve). One gets a good 
approximation over a wide range of non-retarded distances 
0.1 ao < z <C A^. Note how the non-local theory matches 
with the local conductor (solid blue curve) as ao -C z. 

A similar analysis for the CC model [Table Q] yields a 
screening length ao = \J e^D /A-ku involving the bulk dif- 
fusion constant, and the approximate form 



l/a >£/c: r T M(i£, k) 



^k 2 + l/al - k 
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Despite the difference in the reflection amplitudes (I2T1 l23l . 
the short-distance asymptote turns out to be just a factor of 
one half smaller than Eq.(l22l (upper solid gray curve in Fig J2|. 



These results illustrate the dramatic impact of nonlocal elec- 
trodynamics on the van der Waals interaction. We expect them 
to apply qualitatively in materials where the background po- 
larizability f m > 1 provided by bound electrons is relatively 
small. The C3 coefficient in Eq.(fT9b is then nonzero, but weak 
and shifts the short-distance asymptotes of the CL and CC 
curves in Fig|2]below zero. 
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FIG. 2: Atom-surface interaction at zero temperature, from the van 
der Waals into the Casimir-Polder range. We plot the C3 coefficient, 
i.e., the potential U(z) multiplied with z 3 , in units of hQ,a(0). The 
distance z is in units of the screening length ao. The arrow marks 
the (reduced) transition wavelength Aa. Solid (blue) line "local con- 
ductor": Drude metal with dielectric function dl8t , dashed: surface 
covered with a diffusive charge layer (CL model), dotted: continuous 
bulk charge with a hydrodynamic response (CC model). Solid gray 
lines: short-distance asymptotes to the CL and CC models, based on 
the approximate reflection coefficients J21l|23t . 
Parameters: DC conductivity a = 3.6 • 10 17 s -1 (typical for Al), 
dielectric constant too = L electron scattering times r = t s — 
1CP 15 s, diffusion constants D = D s — 5 • 10 3 cm 2 /s, atomic reso- 
nance wavelength 27tAa = 628 nm (Q/2ir = 477 THz). For these 
parameters, the screening length is ao ~ 0.3 nm. The van der Waals 
interaction C3 / z 3 with a local Drude conductor gives a normalized 
value —0.118 from Eq.d20t. The CC/CL models show, at short dis- 
tances, amuch weaker interaction potential 0(1/2) given in Eq.l l22t , 



IV. CASIMIR-POLDER (RETARDED) REGIME 

This regime corresponds to an intermediate range of dis- 
tances, 



A A <C z < A T 



(24) 



where the thermal wavelength A^ = 7.6 at 300 K. This 
range is characterized by the form T ~ 1 /z 4 for both dielec- 
tric and metallic surfaces with a local dielectric response. 

For analytic expansions, we consider T = so that the £- 
integral of Eq.(f3]l still applies. From the exponent 2voz, we 
read off the characteristic frequency tu c = c/(2z) that limits 
the ^-integration range to £ ^ tu c . In the CP regime (1241 1. lu c is 
much smaller than the atomic resonance il, and we can expand 
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the polarizability 



(25) 



The following discussion applies to a free-electron metal at 
high density where = 1 and a ^> f2. 



A. Good conductor 

We then have the small parameter oj c / a to simplify the re- 
flection coefficients. It turns out that the impact of charge 
diffusion is very small in the Casimir-Polder regime: for the 
CL model, we find by inspection that the relevant dimension- 
less ratio is D s / (cz) -C 1. We start with the zeroth order with 
respect to this ratio and expand in powers of 
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Performing the integrations, we get the familiar Casimir- 
Polder potential and next-order corrections 
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where the incomplete Gamma function is T(0,x 
J t~ 1 e~ t dt. The limiting values of third term give ; 
rection -{77/72)5^/2 w -1.34 5 - z" 1 / 2 as uj c t 
and —(8/5)6\/2u! c t w —2.26 8^Juj c t ~ z -1 as lj c t — > oo. 

Note that these approximations correspond to two inequiv- 
alent ways of implementing the perfect-conductor limit. The 
first case could be called "overdamped", with a purely real 
conductivity. The penetration of transverse fields into the 
bulk then occurs by means of diffusion. The correction to the 
Casimir-Polder potential in Eq.(|27T> scales like £ c /z oc 



,-1/2 



where £ c = ^Jc 2 /aui c is the magnetic diffusion length at the 
characteristic frequency oj c . In the second case, the conduc- 
tivity is purely imaginary and, similar to a superconductor, the 
transverse field is screened from the bulk. The correction to 
the Casimir-Polder potential arises from the field penetrating 
a thin layer of the order of the plasma wavelength (also called 
London-Meissner penetration depth) X p = c^J t / (Ana), and 
scales like X v /z. The latter case has been studied, for exam- 
ple, in Ref.[33], Eq.(37), and their result is recovered by the 
two correction terms in Eq.dZTt: 

p w-«y« «4-jft (Aj<< ,«„ )(28) 

Ucp( z ) 3 z A 5 z 

where Ucp(z) is the first term in Eq.(|27]i. This range of 
distances is quite narrow for the parameters of Figj2] where 
cr rj 300 rim. 

Let us now extract the contribution due to the diffusive 
charge layer (CL model). To the first order in the surface dif- 
fusion coefficient D s , the correction to the local model can be 



worked out to be: 

U(z)-U loc (z) ^ D s 5 
U C p(z) ~ 6V^c; 
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In the limiting case lo c t, lo c t s 
a relative correction —(285/88) 
in the opposite limit, — (12/5)a/2t/w 1 



0, the integration gives 

^/2{D s 6/cz) - z- 3 / 2 ; 

2 (D s S/cz) = 



— (6v / 2/5)(c s /c) 2 (A p /z) where the speed of sound c s = 
J D s /t s characterizes the dispersion of longitudinal modes 
in the charge layer. The latter estimate illustrates that the cor- 
rection brought about by the charge layer is negligible in the 
Casimir-Polder regime. A similar conclusion is reached for 
the continuous charge (CC) model; we omit the calculations 
for brevity. The numerical results from Fig J2] illustrate that 
both CL and CC models merge into the local description in 
the Casimir-Polder range z }t \a- 



B. Semiconductor 

Fig[3]shows numerical calculations of the interaction poten- 
tial for a material with a conductivity typical for semiconduc- 
tors. The characteristic length ag for screening is then much 
larger and falls into the Casimir-Polder range of the atom- 
surface potential. The data show that the CC and CL models 
interpolate between the limiting cases of a local Drude con- 
ductor and a non-conducting dielectric (where e(uj) = 
does not diverge at zero frequency). We have taken the rela- 
tively low value = 1.5 to amplify the difference between 
the dielectric and the conductor in the local limit: their differ- 
ence scales with 2/ (e^ + 1). 

In Fig(3j we show the free energy of interaction J-(z) cal- 
culated from the Matsubara sum (|9j. At distances beyond the 
thermal wavelength At, the free energy follows al/z 3 power 
law with a C3 coefficient proportional to T that we discuss 
in the following section. The difference between dielectric 
and Drude conductor arises, for these parameters, from the 
zeroth term in the Matsubara sum. This term is discussed 
in more detail in Sec|V] Indeed, in the other terms, the con- 
ductivity enters only in the ratio Aira/^i = 2Tia/ (Ik-gT). At 
room temperature, this ratio can be neglected compared to the 
background dielectric constant provided the conductivity 
cr <C 4 • 10 13 s _1 . This regime applies to a wide range of 
doped semi-conductors. 

The van der Waals regime for this material is not described 
by Eq.(l20t due to the low conductivity. Ignoring conductiv- 
ity completely, Eq.(fT9b for a local dielectric gives a short- 
range coefficient C3 with a value —1.91 a(0)ksT in the units 
of Fig(3j this corresponds well to the full calculation. We 
have checked that the small difference is actually due to rel- 
atively large deviations from the non-retarded approximation 
that was applied to derive Eq.(fT9b. A similar situation oc- 
curred in Ref . 113411 which discusses the Casimir force between 
two plates separated by a dielectric liquid. 
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FIG. 3: Atom-surface potential through the Casimir-Polder range up 
to the thermal wavelength. We plot the C3 coefficient of the free en- 
ergy, i.e. J r (z)z 3 , but here in units of q(0)£;bT. The distance is 
normalized to the screening length ao « 73 nm. Solid thick curve: 
local dielectric function in Drude form, dashed curve: charge layer 
(CL) model, dotted curve: hydrodynamic (continuous charge, CC) 
model, solid thin curve: non-conducting local dielectric. 
Parameters: background dielectric constant = 1.5, DC conduc- 
tivity a = 10 10 s _1 (comparable to Ge), electron scattering time 
t = t s = 10~ 13 s, diffusion constants D = D s = 4.5cm 2 /s, 
atomic resonance Q./2n — 477 THz (wavelength 2ir A a — 628 nm), 
Temperature T = 300 K (thermal wavelength At = 7.6 /im). 



V. LIFSHITZ (THERMAL) REGIME 

This section deals with the long distance regime At <C z 
where the leading contribution to the atom-surface potential 
is given by the I = term in the Matsubara sum (|9}. The 
other terms are proportional to the exponentially small factor 
cxp(— AttIz/Xt) and can be neglected if the I = term is 
nonzero. A glance at Fig [3] illustrates that the thermal regime 
is already well borne out at z ~ At due to the factor 4V in the 
exponential. 

The static term in the Matsubara sum has been the subject 
of much discussion in the field of dispersion interactions 11351 
3611 . To illustrate this, we give the limiting forms of the free 
energy in the thermal range for an ideal dielectric material 



A T < z : 



a(0)k B T - 1 
Az 3 + 1 



while for a conductor in the same limit 



a{0)k B T 
4^3 



(30) 



(31) 



In fact, the latter result is obtained for any material with a 
nonzero conductivity: as a — > 0, the former (dielectric) result 
is not obtained in a continuous manner fl. This is due to the 
static reflection coefficient 7-tm(0, k) which is equal to 1 for 
any nonzero er, while setting a = from the start for a pure 
dielectric, one gets rxM(0,fc) = (foo — l)/(£oo + 1)- This 
difference between conductor and dielectric is also visible in 
the Casimir-Polder range shown in Fig [3] The discontinuity 



disappears only in the limit T = for the material parameters 
considered here. 

This effect is actually an artefact of the description in terms 
of a local material response (conductivity, dielectric function). 
Using a hydrodynamic model similar to our CC, Pitaevskii 
has shown that the free energy shows a continuous cross-over 
between the limiting cases Eqs. d30l [311 . We show now that 
the same is true for both CC and CL models considered here. 

For the CC model, the first line of Eq.((9) can be written in 
terms of a dimensionless integral (i = 2kz) 



A T < z 



T(z) 



(32) 



a(0)fc B T 



dtt 2 e 



^t 2 + (2z/a a y - t 



^t 2 + (2z/a ) 2 +t 



with the screening length a of Eq. ([]]). We recover Pitaevskii's 
result J3l by calculating ao from the diffusion coefficient 
D rj k-QTr/m* of a non-degenerate electron gas. This leads 
to a,Q 2 = ti£b/(47t) where n is the carrier density in the con- 
ductor and £b the Bjerrum length (i.e., the distance where the 
Coulomb energy between two electrons becomes compara- 
ble to the thermal energy: e 2 /^^^) = k^T). A glance 
at Eq.j32l> tells that the dielectric and metallic values of the 
reflection coefficient are smoothly interpolated as the ratio 
(z/ao) 2 changes from zero to infinity. This is illustrated in 
Fig |4] (dotted line) where the coefficient of the 1/z 3 power 
law is plotted vs z/ao. 

The same qualitative behaviour is found in the CL model 
where the free energy takes the form 



Hz) 



a(0)k B T 
8z^~ 



dtt 2 



: (6 oo -l)t 2 + 6 oo (2z/a ) 2 
(too + l)i 2 + £oo(2z/a ) 2 



(33) 

This is shown in dashed in Fig|4] The parameters chosen here 
are for a very poor conductivity (low carrier density) where 
the screening length ao is large enough to fall into the ther- 
mal range. This applies to dilutely doped semiconductors, or 
to the thermally excited conduction band of an intrinsic semi- 
conductor. 

In Fig [51 we explore under which conditions the atom- 
surface interaction energy is most sensitive to the details of 
the surface charge response. This contour plot shows the ratio 
between the CL and CC results for the interaction energy in 
the thermal range, varying the distance z and the background 
dielectric constant e m . The two models differ maximally in 
the cross over range z ps ao, and for ~ 1. This could 
have been expected from Eqs.d30l [311 because the dielectric 
and metallic limits (the two horizontal lines in Figj4]i are then 
most separated. 



VI. DISCUSSION AND CONCLUSION 

Any material with mobile charges is characterized by a 
screening length ao, and a local current-field relation (Ohm's 
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FIG. 4: Interaction potential in the thermal range: C3 coefficient of 
the free energy, i.e., T(z)z 3 , in units of a(0)kBT. The distance is 
normalized to the screening length ao w 33 fim. Solid thick (blue) 
curve: local conductor (a — 5 ■ 10 4 s _1 , typical for highly purified 
water), dashed curve: charge layer (CL) model with diffusion coef- 
ficient D B = 4.5cm 2 /s, t s = lCF 13 s, dotted curve: continuous 
bulk charge (CC) model (D — D s , r = t s ), solid thin (green) curve: 
ideal dielectric model (eoo = 1.5). The other parameters are as in 
Fig[2] 




(distance z / a ) 

FIG. 5: Contour plot of the ratio between the interaction potentials in 
the CL and CC models. We consider only the thermal range Eqs.( l32l 
[33}. The a>axis gives the ratio (z/ao) 2 , the y-axis the background 
dielectric constant e^. The parameters for surface and bulk diffusion 
are the same (D — D a , r — r s ). 



law) is necessarily limited to scales larger than a . We have 
explored in this paper how screening at and below the sur- 
face influences the long-range van der Waals-Casimir-Polder 
(vdW-CP) interaction between an atom and the body. Our de- 
scription may be termed 'mesoscopic' in the sense that the 
electronic response is collective in nature, but retains traces 



of ballistic carriers via the diffusion coefficient D ~ vi (v 
is a typical carrier velocity and I = vt the scattering mean 
free path). Two models for the electromagnetic response of 
the surface were studied in detail: a continuous charge dis- 
tribution below the surface (CC) within a hydrodynamic ap- 
proximation, and a thin charge layer (CL) with a diffusive re- 
sponse typical for, e.g., localized surface states. Both models 
provide a continuous crossover of the vdW-CP potential be- 
tween two local limiting cases: a pure dielectric and a con- 
ducting medium, as the atom-surface distance goes through 
the range z ~ ao. The two limits can be distinguished from 
the zero-frequency limit of their dielectric functions and give 
different coefficients C3 for the 1/z 3 power laws that prevail 
at very short (van der Waals) and very large (thermal or Lif- 
shitz) distances. Our calculations extend the picture proposed 
in Refs.0,0] to any distances, namely that the nonzero DC 
conductivity a can be neglected if the atom-surface distance z 
is shorter than the screening length ao [Eq.(fl~|i]. 

The differences in the vdW-CP interaction may be used as 
a probe that can identify the type of charge transport in the 
(sub)surface region. The sensitivity of this probe is maximal 
when the screening length ao matches the atom-surface dis- 
tance z, particularly in the retarded range z > A a- TheCC/CL 
difference is particularly large when the dielectric constant 
of bound carriers and background ions is close to 1. The rea- 
son is that the jump — 1 leads to a surface polarization 
charge that responds locally (not by diffusion) and is therefore 
masking the effect of either the CC or the CL. For transitions 
in the optical visible range, these favorable conditions corre- 
spond to z > 1 /im which is indeed the achievable range of 
present experiments. A screening length ao ~ 1 /im occurs in 
a dilute semiconductor, which at room temperature is a non- 
degenerate electron system with ao = \J £ocMbT I (47me 2 ). 
Hence we require a carrier density n w 10 12 cm~ 3 x e^. 
This density can be achieved by dilute doping, as for fused 
silica 0|. One may also work with intrinsic semiconduc- 
tors where the carriers are thermally excited with a density 
n w 10 19 cm" 3 cxp (~E g /2k B T) 0, where E g is the 
gap. Hence, with E g w 1 cV and by varying the temperature, 
one can span a range around 10 12 cm -3 . 

In conclusion, our two charge type models are representa- 
tives of a composite surface of a metal, with a nonlocal elec- 
tromagnetic response due to charge transport in the bulk and 
at the surface. Such composite surfaces are fairly common 
corresponding to either surfaces covered with adsorbates or 
nanostructures, or to disordered surfaces with quantum well 
states. We have studied the Casimir-Polder interaction with 
such surfaces and have shown where this effect can be used as 
a sensitive probe of the surface type and its diffusive proper- 
ties. 
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Appendix A: Surface response 

1. Surface impedances 

The calculation of the electromagnetic Green function pro- 
ceeds by expanding the field created by a point source into 
Fourier components and finding reflection and transmission 
coefficients for each wave vector lq incident on the surface. 
With the wave vector in the xz-plane and the macroscopic 
body in the half-space z < 0, we have lq = ke x — k z e z with 
k z = \Juj 2 /c 2 — k 2 = ivo- We consider separately two prin- 
cipal polarizations. In the TE-polarization, the electric field 
outside the surface is written in the form 



z>0: E(r) = E TE (k)e lkx e 



z ) (Al) 



where e y is the unit vector transverse to the plane of inci- 
dence. One gets the magnetic field from the Faraday equation: 
B x {r) = i( c / w )^z-E'y( r )- The ratio between these two tan- 
gential fields, evaluated at z = + , is the surface impedance 
Zte an d determines the reflection coefficient 



Z 



TE 



Ey(0 + ) 

B x (o+y 



\{cvq/oj)Z' 



TE 



1 



i(cv /uj)Z TE + 1 



(A2) 



In a local model for the body response, one has Zte = 
—iu)/(cv) where the transmitted wave vector is k t = ke x — 
ive z , and the Fresnel formula is recovered. 

In the TM-polarization the electric field vector is in the xz- 
plane, and Eq. dAU becomes 

E(r) = A TM (k)e ikx [(ike z -v e x )e v ° z 

+ r TM (ike z + v e x )e- voz . (A3) 

This gives B y = —i(c/oj)(d z E x — \kE z ) = — i(o;/c)A(l + 
riii) just above the surface. One defines impedance and re- 
flection coefficient from the tangential (x-) component of the 
electric field 



Ztm — 



E x (0 + ) vq-i(uj/c)Z tm 

u 77TTT , r TM = (A4) 

B y (0+) d + i w/c Z T m 



Its local approximation is Ztm = —icv/[uje(uj)] where e(io) 
is the bulk dielectric function. 



2. Solving the reflection problem 

We start to work out the electromagnetic response function 
in the TE-polarization. Within a local description of the bulk 
medium below the layer, one can work with a medium wave 
vector k t in the xz-plane, as defined after Eq. (IA2b . with 



iv = \J (uj/c) 2 e(oj) 



(A5) 



An ansatz similar to Eq. dAU can be written down and aug- 
mented by a longitudinal part 



z<0: E(r) = E TE (k)e ik 



t L (ike x + uie z )e" 12 



(A6) 



where the component V\ of the longitudinal wave vector is as 
yet undetermined. From the Maxwell equations, the tangen- 
tial component E x is continuous, and since it is zero above the 
layer [Eq. dAlH . we find = for the longitudinal ampli- 
tude. The field E z perpendicular to the surface is zero above 
and below the layer, hence the surface charge 7 and the cur- 
rent density J are zero from Eqs.dTO] [TTV The magnetic field 
B x is then continuous as well, and we get the local value for 
the surface impedance from 



Zte = 



Ey(0- 



TE 



B x (0-) i(c/uj)vt T E 
The reflection coefficient ( IA21 > takes the familiar form 

v - v 



rrE 



V + v 



(A7) 



(A8) 



In the TM-polarization, both transverse and longitudinal 
fields in the medium are relevant, as is well known 113 811 . The 
expansion (IA6b becomes 



z < : E(r) = A TM (k)e ikx (ike z - re s )t TM e" 



(ike x + Vl e z )t L e Vl " 



(A9) 



The tangential field E x is continuous and becomes outside the 
layer 



^(0+) = Atm(A:) e ikx (-vt TM + ikt L ) 



(A10) 



Due to the surface current density J x , the field B y has a jump, 
and one gets above the layer: 



47T 

B y (0+) = B y (0~) — —J x 

= — 1 ^TM^TM 



-in 



ikD s (uj)j(AU) 



We need to express 7 and in terms of the transmitted am- 
plitude t.TM ■ use charge conservation (fill and the continuity 
of the z-component of the Ampere-Maxwell equation which 
links the jumps in B y (surface current) and in eE z (surface 
charge). A straightforward calculation yields 



7 



Cr(w)/c.4TMtTM 



(A12) 



uj + \D s {uj)k 2 

and t L = 0. Putting Eqs. JA101 i- (IA12t into Eq.(fA4l). we thus 
find the impedance 

-i(c/cj)u(l + \D s {uj)k 2 /uj) 



Ztm = 



e(u>) + ie 00 L» s (cj)fc 2 /w 



(A13) 



where is the background dielectric function (excluding the 
conduction current). 
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3. Surface conductivity 



Casimir-Polder potential, the estimate (1A181 becomes 



If we model the charge layer as a film of thickness a and 
conductivity cr s {oS), its integrated current density (parallel to 
the layer) has the form: 



J{r,u) = -D s (w)V||7(x,2/) + o s (uj)aE\\ 



(A14) 



where the last term is the conduction current. Including this 
term in the surface response calculations we get reflection co- 
efficients 



v + 4iriaa s (uj)u!/c 



vq + v — At: iaa s {ui)uj / c 2 
with the surface dielectric function being [cf. Eq.dTTIil 

.<r(u) + vaa s ((jj) 



(A15) 
(A16) 



Airi- 



lu + iD s {uj)k 2 



(A17) 



We now identify under which conditions the terms propor- 
tional to <t s (lu) are negligible in these expressions. For an or- 
der of magnitude estimate, we take a s (u>) « cr(u>) and a layer 
thickness a the atomic scale. Therefore we get from Eq. (IA17) 
the condition 



= a\ k 2 - 



-e(w) < 1 



(A18) 



An upper limit can be found easily at imaginary frequencies 
u = i£ where e(i£) = + 4.7rcr/[£(l + £r)] < eoo +W 2 /,!; 2 , 
with Lu p = -J 'Atkj/t the plasma frequency. Using the char- 
acteristic scales k ~ 1/z and £ ~ f2 in the integrals for the 



1 £oo 1 
— + TT + TI < 1 

Z A A A p 



(A19) 



Hence we require the layer thickness a to be much smaller 
than the smallest of the length scales distance z, atomic tran- 
sition wavelength A^ = c/Vl^fe^ and plasma wavelength 
A p = c/ijjp. All these are conditions are well satisfied for a < 
lnm. The surface dielectric functions e' s and e s [Eqs. flA17l ), 
( TTTb l and the ttm reflection coefficients Eq. (IA15b , Tableware 
equivalent. 

As for the tte amplitude dA16l ). we can neglect the surface 
conductivity term provided 



4iraa s (i^ <<; < y 



As an estimate, this is equivalent to 

47rcra 2 £ 



c 2 (l + rO 



< 1 



(A20) 



(A21) 



The maximal value on the left hand side is (a/A p ) 2 , hence 
a must be smaller than the plasma wavelength, as we found 
before in (IA19t . In the limit of large k, the difference vo — v 
in Eq.( IA16b becomes small, and the condition dA20b must be 
replaced by 



kc 2 



(A22) 



Estimating k ~ 1/z, we get a <C z which we also found in 
Eq. (|A19t above. For a < 1 ran, also the tte amplitude is 
therefore not affected by the layer conductivity. 
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